Ion Acoustic Waves in Ultracold Neutral Plasmas 
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We photoionize laser-cooled atoms with a laser beam possessing spatially periodic intensity mod- 
ulations to create ultracold neutral plasmas with controlled density perturbations. Laser-induced 
fluorescence imaging reveals that the density perturbations oscillate in space and time, and the dis- 
persion relation of the oscillations matches that of ion acoustic waves, which are long- wavelength, 
electrostatic, density waves. 



Collective wave phenomena are central to the trans- 
port and thermodynamic properties of plasmas, and the 
presence of a rich spectrum of collective modes is a dis- 
tinctive feature that separates this state of matter from 
a simple collection of charged particles [1]. In ultracold 
neutral plasmas (UNPs) [2|, [3j , which are orders of mag- 
nitude colder than any other neutral plasma and can be 
used to explore the physics of strongly coupled systems 
-0, little work has been done to study collective modes 
jlQj . Here we employ a new technique for creating con- 
trolled density perturbations to excite ion acoustic waves 
(IAWs) in an UNP and measure their dispersion relation. 
This flexible technique for sculpting the density distribu- 
tion will open new areas of plasma dynamics for experi- 
mental study, including the effects of strong coupling on 
dispersion relations [Ill4l4j and non-linear phenomena 
B Oil HH, 03 in the ultracold regime. 

UNPs are formed by photoionizing laser-cooled atoms 
near the ionization threshold. They stretch the bound- 
aries of traditional neutral plasma physics and have ex- 
tremely clean and controllable initial conditions that 
make them ideal for studying phenomena seen in more 
complex systems, such as plasma expansion and equili- 
bration in high-energy-density laser-matter interactions 
[5| and quark-gluon plasmas j6j. UNPs have shown fas- 
cinating dynamics, such as kinetic energy oscillations 
that directly reflect the strong coupling of ions [5|, LLZ[ • 
Strong coupling arises when particle interaction energies 
exceed the kinetic energy |4J]. It is important in many 
fields of phy sics spanning classical to quantum behav- 
ior [3, y, LL8|, LL9| and gives rise to phase transitions and 
the establishment of spatial correlations of particles [4] . 
These studies complement experiments probing strong 
coupling in dusty plasm as 1181 ] and non-neutral plasmas 
of pure ions or electrons [l9j . 

Previous experimental studies of collective modes in 
UNPs were limited to excitations of Langmuir (electron 
density) oscillations with radio frequency electric fields 
[3, [8j which did not determine a dispersion relation and 
were relatively insensitive to dynamics of the strongly 
coupled ions. A high-frequency electron drift instability 
was observed in an UNP in the presence of crossed elec- 
tric and magnetic fields [9[. Spherically symmetric ion 
density modulations were shown to excite IAWs in nu- 
merical simulations of UNPs [10j . Here we excite IAWs 
through direct imprinting of ion density modulations dur- 



ing plasma formation and image them in situ with time 
resolved laser-induced fluorescence [20[ . 

Low frequency electrostatic, or longitudinal ion density 
waves are one of the most fundamental oscillations in a 
plasma along with Langmuir oscillations [2l|. A hydro- 
dynamic plasma description, assuming slow enough ion 
motion for electrons to remain isothermal and an infinite 
homogeneous medium, predicts the dispersion relation 
for frequency uo and wavevector k [l| 
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where M is the ion mass, T e is electron temperature, and 
Xd = \f€oT e /n e e 2 is the Debye screening length, for elec- 
tron density n e and proton charge e. We have neglected 
an ion pressure term because the ion temperature sat- 
isfies Ti <C T e in UNPs. In the long wavelength limit, 
which is the focus of this study, th is mode takes the form 
of an I AW with uj = k^/kBT e /M, in which ions provide 
the inertia and electrons provide the restoring pressure. 

IAWs are highly Landau damped unless Ti <C T e [l|, 
however, they have been studied in many types of high- 
temperature laboratory plasmas pjj Us, [221 . Closely 
related to this work, acoustic waves of highly charged 
dust particles in dusty plasmas have been studied ex- 
perimentally [23|, l24[ and theoretically |lll4l4j| because of 
the possibility of observing the effects of strong coupling 
on the dispersion relation, but to date these effects have 
been masked by damping due to collisions with back- 
ground neutral gas pjj [l3|, [25| . Beyond fundamental in- 
terest, IAWs are invoked to explain wave characteristics 
observed in Earth's ionosphere [26| and transport in the 
solar wind, corona, and chromosphere [27[. 

UNPs are created through photoionization of laser- 
cooled strontium atoms from a magneto-optical trap 
(MOT) [3]. The MOT operates on the 88 Sr ^o^Pi 
resonance at 461 nm [30], trapping ~3x 10 8 atoms at 
~ 10 mK with a spherical Gaussian density distribution. 
Photoionization is a two-photon process: the first from 
a pulse amplified laser beam operating on the Sr trap- 
ping transition, and the second from a 10 ns Nd:YAG 
pumped dye laser tunable around 412 nm [3]. This pro- 
cess ionizes ~50% of the atoms, and the plasma inherits 
its density distribution from the neutral atoms. Effects 
of un-ionized atoms on the plasma are not observed be- 
cause of the fast time-scales of the experiment and small 
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FIG. 1: Excitation, imaging and modeling of density 
perturbations in an UNP. a, A periodic transmission 
mask modulates the ionization beam intensity to excite 
IAWs. An adjustable delay later, a sheet of laser light 
near resonant with the principle transition in the ions, 
illuminates a central slice of plasma, and laser-induced 
fluorescence along a perpendicular direction is captured 
by an imaging system to form a two-dimensional false 
color plot of the ion density distribution, b, ID slice 
through the density profile for T e (0) = 105 K and 1 mm 
mask period, 750 ns after ionization. Deviations from 
the fit Gaussian represent the IAW density modulation, 
Sn. c, Evolution of Sn for the same initial conditions as 
figure [Ud. Time since ionization is indicated on the 
right, and Sn has been scaled by instantaneous peak 
density no(t) and offset for clarity. Solid red lines are 
fits to equation (j4j), and the dotted black lines follow 
the wave nodes and antinodes with time. 



neutral-ion collision cross-sections. 

Resulting electron temperatures are determined by the 
excess energy of the ionizing photons above threshold, 
and are adjustable from 1-1000 K. Ion temperatures of 
approximately 1 K are set by disorder-induced heating 
[171 . l28l |29] during the thermalization of the ions, and 
result in strongly-coupled ions in the liquid-like regime 0, 
0|. Density distributions are spherical Gaussians, n(r) = 
noexp(— r 2 /2ao), with no ~ 10 15 cm -3 and <Jo ~ 1.5 mm, 
yielding average A^ from 3-30 jam. In our experiments, 
no magnetic field is applied. 



IAWs were excited by passing the 412 nm ionizing 
beam, after its first pass through the plasma, through 
a periodic transmission mask and retroreflecting it back 
onto the plasma. This creates a ~ 10% plasma density 
modulation with wavelength set by the period (Ao) of 
the mask (Fig. []Ji). The mask pattern is translated to 
align a density minimum to the center of the plasma. 
Small higher-harmonic IAWs, arising from the square- 
wave nature of the mask, were observed for longer period 
gratings, but no effect on the fundamental wave was de- 
tected. 

For a diagnostic, ions are optically excited on the pri- 
mary Sr + transition, 2 Si/ 2 — 1 Pi/2, with a tuneable, nar- 
rowband (~ 5 MHz) laser at A = 422 nm, propagating 
approximately perpendicular to the ionizing laser [3|, [20j . 
The 422 nm beam is masked with a 1 mm slit so it only 
illuminates a central slice of the plasma, and resulting 
laser-induced fluorescence emitted close to perpendicular 
to the plane of the ionizing and 422 nm beams is im- 
aged onto an intensified CCD camera with a resolution 
of 13 jam. 

Data from ~ 50 repetitions of the experiment is 
summed to form a single image, F(x,y,u), that has a 
frequency (y) dependence reflecting the natural linewidth 
and Doppler-broadening of the transition [3, l20[ . 40 im- 
ages are recorded at evenly-spaced frequencies for the 
422 nm laser spanning the full spectral width of the sig- 
nal. These are summed to obtain a signal proportional 
to the density of the plasma in the illuminated plane 
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Density averaging along the imaging axis (z) is small be- 
cause fluorescence is only excited in a sheet of plasma. 
Absolute determination of density is obtained by cali- 
brating fluorescence signals against absorption images of 
the plasma [3, 20]. 

Figure[T]3 shows a central ID slice through density data 
along the modulation direction. The difference between a 
fit Gaussian distribution and the data yields the density 
perturbation Sn as shown in Fig. [1J. Note the oscillation 
of the wave in time and the changing wavelength of the 
excitation. 

The perturbation has a negligible effect on the global 
dynamics of the plasma, which is dominated by expan- 
sion into the surrounding vacuum. For an initial Gaus- 
sian density distribution and conditions such as used here 
in which inelastic electron- ion collisions are negligible, 
the expansion is self-similar and the characteristic plasma 
size changes according to J3[ 
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where r exp « y/Mo^JkBTJfi) is the characteristic 
plasma expansion time, which ranges from 10-30 jis in 
this study. Associated with expansion is an adiabatic 
cooling of electrons according to T e (t) = T e (0)/(1 + 



t 2 / r exp)i which shows that expansion is driven by a trans- 
fer of electron thermal energy to the kinetic energy of 
ion expansion. In direct measurements of plasma size 
and spectral measurement of ion expansion velocity [20| , 
we saw no deviation in expansion dynamics due to the 
presence of density modulations. 

To describe wave data such as Fig. QJ, we scale the 
amplitude to the peak density at t = of the Gaussian 
fit (Fig. [lb]) and assume a damped standing wave with a 
Gaussian envelope, 
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This expression fits the modulation at a single time, t, 
and yields the instantaneous amplitude A(t), envelope 
size a env (t) and wavevector k(t). The hydrodynamic de- 
scription for an infinite homogeneous medium used to 
derive the IAW dispersion relation, Eq. [TJ allows for pla- 
nar standing- wave solutions, but finite size, plasma ex- 
pansion, and density inhomogeneity introduce additional 
factors that are not small here and preclude an analytic 
solution, so this model is only phenomenological. 
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FIG. 2: Evolution of the IAW wavelength and plasma 
size, normalized to initial values, for T e (0) = 70 K and 
<Jo = 1.45 mm. Both quantities follow a universal curve 
indicating the wave vector expands with the plasma. 
For the solid line, r exp has been set to its theoretical 
value. 



To examine the effect of the expanding plasma on the 
wave we compare the change in the wavelength with 
the changing size of the plasma by plotting, in Fig. [2j 
X(t)/Xo = ko/k(t) and <t(£)/<to, where all quantities have 
been normalized to the values at t = 0. Initial wave- 
lengths match the period of the mask used. All data 
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The agreement for the wavelength indicates the wave is 
pinned to the expanding density distribution. 

To extract the frequency of the wave, u(t), we fit the 
amplitude variation in time, accounting for the evolution 
of the phase, as 



(/' 



A(t) = A e~ rt cos [</>(£)] = A e~ rt cos ( / tj(t')dt' 

'(5) 

We assume the form of the dispersion for an infinite, ho- 
mogeneous medium from Eq. [TJ the observed variation 
of wave- vector k(i), and electron temperature evolution 
T e (t) predicted for a self similar expansion [3[ to obtain 



Lj(t) = k(t)y / k B T e (t)/M = UJo 
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Initial amplitude Aq and frequency Co>o , and damping rate 
T are allowed to vary in the fits, which match the data 
very well (Fig. [3]). Decreasing frequency with time was 
observed in simulations of spherical IAWs [lOj . 
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FIG. 3: Evolution of wave amplitude and fits to obtain 
cj(t). The main plot shows the instantaneous amplitude 
for various mask periods with T e (0) = 70 K. The inset 
similarly compares data for different initial electron 
temperatures with Ao = 0.50 mm. Lines are fits to 
equation ([5]) in which r exp has been fixed to the 
theoretical value. 



For a range of initial electron temperatures and mask 
periods, we extract uo and fco an d calculate the disper- 
sion of the excitations, as shown in Fig. [U The excel- 
lent agreement with theory, Eq. [lj confirms that these 
excitations are IAWs. The planar standing- wave model 
captures the dominant behaviour of the wave, and to a 
high accuracy there is no deviation from the standard 
dispersion relation in spite of the plasma's finite size, ex- 
pansion, and inhomogeneous density. 

Following on this initial study of ion density excita- 
tions in an UNP, there are many topics to explore. The 
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FIG. 4: Dispersion relation of IAWs for different initial 
electron temperatures. Lines are from the theoretical 
dispersion relation, Eq. [TJ with no fit parameters. 



observed waves show damping times on the order a few 
oscillation periods, which is faster than predicted for Lan- 
dau damping [l[ . The nature of the boundary conditions 
and effects of density inhomogeneity and plasma expan- 
sion on IAWs should be explored, and may be important 
for understanding damping. The mask period is cur- 
rently limited by diffraction effects, but this can be over- 
come by modifying the optical configuration. For mask 
periods approximately five time smaller, we can probe 
beyond the acoustic region of the dispersion relation. In 
this regime, the effects of strong coupling on dispersion 
are also predicted to be important, and there are many 
theoretical predictions that have not been tested [Ill4l4| . 
The wave can be studied in velocity space with resolved 
fluorescence spectroscopy [20], and different initial den- 
sity distributions can be designed to investigate solitons 
[15j, instabilities [16j, asymmetric excitations, and shock 
waves [31 KLOj. 
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